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The fundamental role of quantized vibrations in coherent light harvesting by 
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The influence of fast vibrations on energy transfer and conversion in natural molecular aggregates is an issue 
of central interest. This article shows the important role of high-energy quantized vibrations and their non- 
equilibrium dynamics for energy transfer in photosynthetic systems with highly localized excitonic states. We 
consider the cryptophyte antennae protein phycoerythrin 545 and show that coupling to quantized vibrations 
which are quasi-resonant with excitonic transitions is fundamental for biological function as it generates non- 
cascaded transport with rapid and wider spatial distribution of excitation energy. Our work also indicates 
that the non-equilibrium dynamics of such vibrations can manifest itself in ultrafast beating of both excitonic 
populations and coherences at room temperature, with time scales in agreement with those reported in 
experiments. Moreover, we show that mechanisms supporting coherent excitonic dynamics assist coupling 
to selected modes that channel energy to preferential sites in the complex. We therefore argue that, in the 
presence of strong coupling between electronic excitations and quantized vibrations, a concrete and important 
advantage of quantum coherent dynamics is precisely to tune resonances that promote fast and effective energy 
distribution. 
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I. INTRODUCTION 

The functionality of natural multichromophoric com- 
plexes depends fully on the strength and structure of in- 
teractions between electronic states and vibrational de- 
grees of freedom associated either to intramolecular vi- 
brations or to intermolecular interactions with the pro- 
tein and the surrounding solvent environmeniP^. In 
these systems, the phonon environment influencing elec- 
tronic excitation dynamics can be broadly divided into 
two distinct components^. Firstly, a continuous distribu- 
tion of modes with a characteristic cutoff frequency below 
or comparable to the thermal energy scale, ksT. These 
modes, arising from the solvent and low-energy protein 
vibrations, induce thermal fluctuations of onsite ener- 
gies which in turn drive transitions between energetically 
close excitonic states. Secondly, spectroscopy studies^^^ 
also reveal the active participation of specific vibra- 
tional modes during excitation dynamics, often described 
as peaks in the spectral density characterising exciton- 
phonon interactions. Vibrations with lower frequencies 
may arise from inter-molecular coupling between pig- 
ments and the local protein environment^^ while peaks 
at high energies originate mostly from intramolecular vi- 
brations that give rise to vibronic transitions^^EH^ \Ye 
denote these high-energy vibrational modes as quan- 
tized vibrations to indicate that their quantum me- 
chanical features can be significant even at ambient 
temperatures^. In contrast, vibrations with frequencies 
below ksT will be thermally activated. Indeed, the in- 
fluence of quantized vibrations in energy transport in a 
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variety of nat ural mo lecular systems is currently of cen- 
tral intereslP«SE3IIl. 

In this work we discuss the important role of quan- 
tized vibrations and their non-equilibrium dynamics for 
electronic excitations in the cryptophyte antennae pro- 
tein phycoerythrin 545 (PE545)^^^ where exciton co- 
herence beating has been observed^. Chromophores in 
cryptophyte light-harvesting antenna proteins have large 
energy gaps and are separated by an average distance of 
20 A^^, which is double the average pigment distance in 
proteins of plants^, hence the majority of electronic cou- 
plings in PE545 are small (<30cm~^p^. How can this 
antenna protein function efficiently with such structure? 
We show that a fundamental mechanism for the biologi- 
cal function of this complex and for supporting coherence 
oscillations at long-times is the coupling between excitons 
and vibrations quasi-resonant with excitonic energy dif- 
ferences. The importance of these interactions in similar 
light-harvesting systems has recently been highlighted by 
ultrafast spectroscopy experiments^^ "^^. 

Vibrations r eporte d in spectroscopy of PE545 have 
high-frequencies^^ ^^ with the largest couplings for modes 
between 514 cm~^ and 1450 cm~^. Given that these 
frequencies match closely energy gaps between exciton 
states, transport in PE545 is fundamentally activated 
by these vibrations giving rise to the transfer pathways 
and transfer times observed in experiments^. This result 
can be anticipated using Forster theor}^^^ but is derived 
here using a non-Markovian polaron-represe ntatio n mas- 
ter equation for multichromophoric system^^^. This 
more general framework allows us to investigate the evo- 
lution of excitonic coherences and to show that non- 
equilibrium dynamics associated to such vibrations can 
contribute to the long-time oscillations of the coherences 
probed in recent ultrafast experiments^^'^. Based on 
this analysis, we are able to elucidate two key features of 



excitation transfer assisted by vibrations in crypt ophytes. 
First, excitation dynamics follows a non-sequential order 
bridging large energy gaps and distances and therefore 
generating faster and wider spatial distribution of excita- 
tion energy. Second, vibrations that modulate the evolu- 
tion of excitonic coherences also promote energy transfer 
to preferential sites in the complex. Both of these process 
support biological function of the antennae by assuring 
that energy transfer to other photosynthetic complexes 
happens in the relevant biological time scale. On the ba- 
sis of these results we argue that in complexes with highly 
localized excitons a fundamental role for quantum coher- 
ent contributions to dynamics is precisely to ensure that 
resonances with selected vibrations are tuned to guaran- 
tee faster and effective distribution of excitation energy. 
The interactions between excitons and quantized vibra- 
tions in close resonane with exciton gaps are therefore an 
important design principle for optimal light-harvesting 
antennae^. 



II. MARINE CRYPTOPHYTE ANTENNAE PROTEINS 

Cryptophyte photosynthesis relies on phycobiliproteins 
such as PE545 and Phycocyanin 645 (PC645) to en- 
sure photon absorption in the visible part of the spec- 
trum, with an optimal performance in the blue (545 



PE545 given by {h = 1): 



nm) and green (645 nm) region j * I Excitation en- 
ergy absorbed by cryptophyte antennae proteins is trans- 
ferred towards photosystems II and II (PSI and PSII)^^. 
Crystallograph}B^^ reveals that PE545 contains eight 
light- absorbing bilin molecules covalently bound to a four 
subunit protein scaffold labelled A, B, C and D (see 
Figure la). The lowest-energy chromophores are 15,16- 
dihydrobiliverdin chromophores and are linked to sepa- 
rate protein subunits A and B (DBVa and DBVb). The 
remaining six chromophores are phycoerythrobilin (PEB) 
bilins, two of which are doubly covalently bound to the 
protein environment while the remaining four are singly 
bound. The PEB bilins are arranged such that pigments 
PEB50/6IC5 PEB82C and PEB158C occupy the subunit C, 
while chromophores PEB50/6ID7 PEB82D and PEB158D 
are attached to subunit D. One particular feature of this 
complex is that it contains a deep, water- filled slot ex- 
posing the pigments to a large contact with water. The 
combined pigment-protein structure is shown in Figure 
la, while the relative positions of the isolated pigments 



is illustrated in Figure lb 
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+ V Knn(cr+Cr^ ^^n^m)^ (1) 



{m^n) 



with (j^ being the operator that creates an excitation on 
site m i.e. cr+|0) = \m). A salient characteristic of this 
system is that site energy differences \em — e^l are very 
large in comparison to inter-site electronic interaction 
Vmm for instance, the largest coupling is V48 = 92cm~^ 
between the central pigments PEB50C and PEB50D while 
their energy difference is about 1040 cm~^ (remaining 
values of Vmn can be found in the supplemental infor- 
mation). As a result, each excitonic eigenstate, which 
we denote |e^), is highly localized to a particular site 
m. The table in Figure Ic presents the site energies of 



the eight chromophores in ascending order as well as the 
subunits to which each of them is bound. Notice the 
interesting relationship between the spatial arrangement 
of the pigments and their energies: the three highest en- 
ergy chromophores are located in subunit D, the next 
three highest energy chromophores lie in subunit C and 
finally the two lowest chromophores reside in subunits 
B and A. Therefore, the hierarchy of pigment energies 
imposes also a hierarchy of energies among protein sub- 
units. In section [Vl we illustrate how one of the functions 
of quantized vibrations is precisely to support efficient 
energy transfer between protein subunits. 



B. System-Environment interaction 

Optimal energy transfer in light harvesting systems de- 
pends crucially on the transitions and correlations that 
can arise from interaction between the electronic sys- 
tem and their environment. The later is usually de- 
scribed as a bath of harmonic modes of frequencies 
cjk, with boson creation operators h\^ and Hamiltonian 

Hb = Xlk'^k^k^k- The linear interaction of this har- 
monic bath and the electronic system is given by 



Hi 



^<^m<^n 
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'5k,m^k), 



(2) 



A. The Electronic System 



where ^k,m is the coupling of site m to mode k. The 
strength and structure of these interactions is character- 
ized by the spectral density, J{uj) = ^j^ \g\^\'^S{uj — cjk) 
which has two contributions: coupling to low- frequency 
continuous modes denoted as Jcm{^) and coupling to 
quantized vibrations denoted Jqm{^) such that J{uj) = 

JCM (^) + JqM (^) . ^^ 

Guided by linear spectra studies on PE545^^^, in this 
work we consider that Jcm{^) is given by a super-Ohmic 
spectral density of the form 



Quantum chemical calculations^^^ as well as 
spectroscop}^^^*^ have allowed the refinement of a 
parametrized Hamiltonian for the electronic structure of 



JcmH - ^e 



2^1 



(3) 



Here the parameters (in units of cm ^) for the contin- 
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(a) Model of PE545 showing the four protein subunits 
as colored ribbons 



(b) Structure of the isolated pigment molecules in 
PE545 
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(c) Table showing bare excitation energies for each pigment and the protein subunit to which it is bounded. 
These values have been taken from Ref.^^ and they correspond to a model that best fit experimental results. 

FIG. 1. Model of an isolated PE545 from cryptomonad Rhodomonas CS24 based on the X-Ray structure at 0.97A resolution 
(Protein Data Bank ID code IXGOp-^ 



uous modes are si = 3.46 x IQ-"^, 52 = 2.02 x IQ-"^, 
Qi = 1.45 X 10-^ and ^2 = 4.34 x 10"^. We have chosen 
a super- Ohmic form for the continuous component of the 
spectral density taking into account that the frequency 
dependence of this function has been shown to appropri- 
ately reproduce vibrational sidebands^^. Moreover, the 
theoretical framework that we will use fo r describing dy- 
namics is more accurate in this case^^^. 

SpectroscopySni q\^q reveals the influence of 14 vibra- 
tional components in the excitation dynamics on PE545, 
whose frequencies ujqm and Huang- Rhys factors sqm are 
shown in Table |Tj Notice that the majority of the vi- 
brations satisfy the condition Hujqm ^ ksT ^ 200cm~^ 
falling in the regime of high-energy. The largest cou- 
plings (given by ^ = ujqm^/sqm) are for modes ujqm = 
1111 cm~^ and ujqm = 1450 cm~^. In this work we 
include the effect of these modes via a spectral den- 
sity exhibiting sharp components. When combined with 
the non-Mar kovian polaron treatment (see III), this pro- 
vides a starting point to investigate the effects of non- 
equilibrium dynamics associated to these vibrations. In 
practice, these modes could be broadened by interactions 
with the surrounding environment. Therefore, in this 



work we assume broadened vibrational modes, each with 
a Lorentzian line shape of broadening iP'^ . The contri- 
bution of each of these modes to the spectral density is 
given by: 



Jqm{^) 



2sqm^qm 



j% 



(Cj2 



'^QM 



)2 + b\ 



(4) 



A flrst insight into the relevance of the different compo- 
nents of J{uj) in excitation dynamics can be gained by 
understanding the relative position of excitonic transi- 
tions to the spectral density. Figure [2a| shows the forms 
of Jcm{^) and Jqm{^) as well as superimposed bars in- 
dicating the excitonic energy differences of the the elec- 
tronic Hamiltonian (Eq. fl])). One striking feature in this 
flgure is that the energies of the vibrations closely match 
a large number of the excitonic transitions. For exam- 
ple, the localized mode at uqm = 1111 cm~^ lies close 
to resonance with the transition between eigenstates \es) 
and 1 64) and also to the transition between eigenstates 
1 67) and 1 62). To understand the physical signiflcance of 
these resonances, consider a scenario where a single lo- 
calized mode is exactly on resonance with an excitonic 



sqm 



Sqm 



Sqm 

^QM 



0.0013 0.0072 0.045 0.0578 0.045 
207 244 312 372 438 



0.0924 0.0761 0.0578 0.0313 0.0578 
514 718 813 938 1111 



0.1013 0.0265 0.0072 0.0113 
1450 1520 1790 2090 



TABLE I. Vibrational parameters for PE545 taken from 
Ref.^^. Frequencies in units of cm~^. 



transition. In this case, the electronic system is able to 
emit a phonon into the mode, preserving energy dur- 
ing the transition. Therefore, quasi-resonant vibrations 
can activate otherwise energetically unfavorable excitonic 
transitions, bridging large energy gaps and, in this case, 
also large distances. 

Environmental Correlations: The mechanisms 
supporting beating of excitonic coherences in light- 
harvesting systems are still under and they can indeed 
be system-dependent. In the case of PE545, it has been 
argued that non-Markovian dynamics arising from slow 
relaxing modes or solvent-induced correlated fluctuations 
of on-site energies can help supporting long-lasting oscil- 
lations of excitonic coherences'^. So far however none of 
these hypothesis has been experimentally probed. In- 
deed, atomistic descriptions have challenged the idea 
of spatially correlated fluctuations in other molecular 
complexe^L However, solvent-induced correlated fluc- 
tuations can be expected to play a more significant role 
in systems where chromophores are in close proximity 
to water. An example is hydrated DNA where hole 
transport is affected by correlated fluctuations that arise 
from the nearby solvent^^. Interestingly, the structure of 
PE545 is such that a deep pocket of bound water occupies 
space between the monomers, thus allowing close contact 
of the pigments with the water molecules. Atomistic cal- 
culations to understand whether interaction with near- 
by solvent can induce correlated fluctuations in PE545 
are under development^. Therefore, in this work we 
compute linear spectra and dynamics in PE545 with and 
without solvent-induced correlated fluctuations. 

Spatially correlated fluctuations are introduced via 
position dependent coupling to common propagating 
modes, i.e. ^k,m = gi^e^^-^"^ wh ere V m denotes the po- 
sition of the 771 'th chromophore^ffiS Generally, this 
behaviour is limited to wavelengths on which the en- 
vironment appears homogeneous, restricting the modes 
to the low energy sector. Therefore, we assume that 
only the continuous component of the spectral den- 
sity induces correlated fluctuations. For an isotropic 
three-dimensional bath, the spectral density function 
of fluctuations between sites m and n can be written 
as J--(a;) = J^-{uj) + J^- (c^), where J^-{uj) = 

Jcm{^) sinc(u;r^n/'^p/i) and Jqm{uj) = JQM{^)5mn' 
Here Vmn corresponds to the distance between the m'th 



and n'th chromophores, while v^h = 1000 ms~^ is com- 
parable to the speed of sound in water. If instead Vph 
is set to zero, we recover the case of independent har- 
monic baths. For a detailed justification of this model of 
solvent- induced correlations see Nalbach et olW^. 



III. DYNAMICS AND SPECTRA IN THE 
POLARON-REPRESENTATION 

Exact treatment of the exciton-phonon interactions 
across the whole rang e of energ y and couplings is not 
a straightforward task^^^^^^^^H^. Importantly, in the 
regime of strong system-environment coupling, Forster 
theory^ accounts well for the influence of high-energy vi- 
brations in transfer rates. However, Foster theory fails to 
capture the non-equilibrium dynamics of such vibrations 
due to interactions with the electronic components. A 
formal treatment of quantized vibrations consider them 
as part of what is denoted as system} such that coher- 
ent quantum interactions between electronic and vibra- 
tional modes are accounted for. We illustrate such an 
approach with an example in section VC but the full 
quantum treatment of PE545 with 14 vibrational modes 
per pigment is beyond the scope of this paper. In- 
stead, here we use a non-Markovian many-site polaron 
master equation^^ ^^, which allows consideration of non- 
equilibrium dynamics of high-energy quantized vibra- 
tions, that is, the fact that these modes take some time 
to thermalize. As we will see here, these non-equilibrium 
dynamics induce modulation of excitonic population and 
coherences. The treatment is particularly accurate for 
systems where excitonic states are highly localized and 
hence suitable for describing dynamics and spectroscopy 
of PE545. Full details of the polaron master equation for 
multichromophoric systems are given in Refs.^^^^. No- 
tice that this formalism allows reconstruction not o nly of 
the site-population dynamics as originally proposecP^I^ 
but also of the evolution of excitonic coherences^. For 
self-consistency of the paper we summarize the frame- 
work here. 

The polaron formalism describes a dynamics in a 
framework of displaced oscillators where the original 
Hamiltonian H = H^i + Hb + Hj is transformed via a 
small polaron-transformation i.e. H = e^He~^ with S = 

Em^m^mEk(^k,m^k " ^k,m^k) and /lk,m = ^k,mM. 

In this shifted-oscillators framework, the effective elec- 
tronic system is described by 



/ ^ ^m^rn^rn 



E 

{m,n) 



^mnPmn\^rn^n ^ ^n ^mJ' 



(5) 



Here Cm and VmnPmn are shifted onsite energies and 
renormalized electronic interactions respectively. The 
eigenstates of Hei are labelled |a^) and are denoted as 
renormalized exciton states from now on. Interaction 
with the bath of shifted oscillators induce fluctuations of 
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(b) Schematic of prototype pairs in PE545 



FIG. 2. (a) Main panel shows the contribution of quantized vibrations Jqm{oj)^ while the inset depicts the continuous spectral 
density Jcm{(jo). Bars indicate excitonic transitions, (b) Some pigment pairs have energy detunings A^ much smaller than 
any of the cjqm, for instance the pair DVB a and DVBb (sites 1 and 2) for which A^ = 35 cm~^, thus quantized vibrations 
have a small effect on this excitonic transition. In contrast, many pairs formed of pigments located on different protein subunits 
have AE in close resonance with a vibrational mode, example of this is the central dimmer PEBsoc and PEB50D (sites 4 and 
8) whose energy difference is 1040cm~^. In this later case, electronic interaction Vmn is the smallest energy scale. 



electronic couplings given by 



Hi 



E 

{m,n) 



Vmn{B, 



mn^rn^n 



Bi 



mn^n ^m) 



(6) 



Bath-induced renormalisation factors are defined as 
l^mn = (Bmn) while displaced bath operators become 

B — B — B with B — pEk(^^k,rr.ri&l,-5/lk,rr^ri^k) 

When a common bath is assumed, the factor (^/ik,mn = 
^k,m — ^k,n depends on the difference in electron-phonon 
couplings of sites m and n. For a harmonic oscillator 
bath in thermal equilibrium, the renormalisation factors 
evaluate to /3^n = e"^ ^kCoth(/3u;k/2)|5ak,^.P^ 

Perturbation theory with respect to Hj^ leads to a non- 
Markovian polaron master equation 



dt 



~p{t) = n{t)~p(t) ^x(typ{{)), 



(7) 



The homogeneous term 1Z{t)p{t) describes the relaxation 
of the electronic degrees of freedom due to the interaction 
with a thermalised environment. The inhomogeneous 
term X{t)p{0) is only non-zero when the initial environ- 
mental state differs from the thermal state in the polaron 
frame. Therefore, this term allows us to capture the non- 
equilibrium dynamics of all vibrational components and 
in particular the fact that high-energy quantized vibra- 
tions are taken out of thermal equilibrium due to inter- 
action with the electronic componenlP^l^. Explicit ex- 
pressions for the homogeneous and inhomogeneous terms 
are presented in the supplemental information. 

Original frame observables The master equation 
presented in Eq. (JTl) describes the evolution of system 
observables within the polaron frame {p{t)). However, 



we are ultimately interested in the dynamics as seen in 
the original untransformed 'lab' frame {p{t)) since this 
will be the appropriate operator to compute expected 
values of physical observables. A general prescription 
for the transformation of system observables from the 
polaron to the lab frame was presented in Kolli et al.*^. 
These transformations can be summarised as follows: site 
populations (cr^cr^) remain unaffected upon transforma- 
tion to the lab frame as the operator cr^^cr:^ commutes 
with the polaron transformation S. Expected values of 
off-diagonal operators in the site basis are given by an 
elaborated transformation, which up to first order in Hj 
reads: 

i^m^n) = tr{cr+cr-p(0)} + /3^ntr{cr+ a" [p(t) - p(0)]} 
+tr{a+a-5i(t)p(t)} + tr{a+a-Ti(t)p(0)}. 

Expression for super-operators Si{t)p{t) and Ti{t)p{0) 
are presented in the supplemental information. 

Linear Spectra: Here we provide a derivation of the 
linear spectra based on the polaron-representation of the 
dynamics. Our derivation of the spectra appeals to the 
Markov and secular approximations in the polaron frame, 
formulation which may still capture some non-Markovian 
effects in the untransformed frame (see details in Kolli et 
al?^). The absorption spectra of a multichromophoric 
system is defined in terms of Fourier transform of the 
dipole-dipole correlation function^^: 



A{ijo) =Re dt e'^^ Cd-dif) 



(8) 



The dipole-dipole correlation function is defined as the 
two-time correlation function of the transition dipole op- 



erator 



Cd-d{t) = Tvs+B{Kt)fi{0)pB\0){0\} 



(9) 



where /i(t) = e~ fie is the system transition dipole 
moment in the Heisemberg picture. The system tran- 
sition dipole operator fi is given by the sum of the 
dipole operators for the individual chromophores ji = 
^rn^l^'rn\^) {'^\^h.c.) . The dipole-dipole correlation func- 
tion after transformation to the polaron frame reads 



Cd-d{t) = Trs+B{A(t)M(0)pB|0){0|} 



(10) 



where the polaron frame dipole operator is given by: 
M = Xin (^(-^k,m)Mm|0)(m| + /i.e.) 

m k 
Z! n (^(-/ik,m)MmWma|0)(a| + h.c)j (11) 



a.^m k 



where I^(/ik,m) = e^- 



k,mOi^-All^^^Ok 



is the bath displace- 



ment operator of mode k due to interaction with site m 
and Uma = {Tn\a) is the amplitude of site m on renor- 
malized eigenstate \a). 

Under the Markovian and secular approximations 
within the polaron frame, the dipole-dipole correlation 
function can be written as: 



Cd-d{t) = ^Ca{t)Ga{t)eM^^at) 



(12) 



For each renormalized exciton \a) there are two contribu- 
tions to the dipole-dipole correlation function: one aris- 
ing from the displacement of the vibrational environment, 
and a component due to dephasing of the electronic de- 
grees of freedom within the polaron frame. The bath 
correlation function in the polaron frame Ca {t) is defined 
as 



Ca{t) = ^l^ml^nUmaUnaT^'^B{D{-h\^^n)Dih\^,mit))pB} 

mn 

(13) 
The electronic correlation function is Ga{t) = exp(— r^^t) 
with a dephasing rate for each transition given by 

f =vrr(^) +r(^) +r(^) +r('^) ^ (u) 

and r^*i Markovian polaron rates given in Kolli et 
al.^. Full details of the above derivation is presented 
in Appendix [Aj 

The final form of the absorption spectra is given by 

POO 

A{uj) = V i?e / dt e*("-"°)* Ca{t)Ga{t) 

(15) 
In a similar manner we can derive the following expres- 
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(a) PE545 absorption spectrum at room temperature (294 K) 
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(b) PE545 fluorescence spectrum at room temperature (294 K) 



FIG. 3. Comparison of experimental (solid black) and cal- 
culated (solid red) spectra. Dashed-dotted curve indicates 
calculated spectra in the absence of solvent-induced correla- 
tions. The results presented include average over Gaussian 
static disorder for the site energies with FWHM= 500 cm~^ 



sion for the fiuorescence spectra: 

POO 

F{uj) = Tn{uj^)Re / dt e'^'^"'^")* C:(t)G;(t) 

(16) 

where n{uja) = i^^viP^a) ~ ^)~^ is the thermal occupa- 
tion of eigenstate a. 



IV. SIGNATURES OF QUANTIZED VIBRATIONS IN 
SPECTRA AND DYNAMICS OF PE545 

A. Signatures of quantized vibrations in the linear spectra 

Experimental signatures of specific vibrations partic- 
ipating in excitation dynamics in light-harvesting sys- 
tems can be seen in photon-echo spectrocopy^, reso- 
nance Raman studies^, fluo rescence narrowing^ and 
hole-burning experiment s^Eoj. Here we show that the the- 
oretical framework considered reproduces known room- 
temperature spectra and we point out signatures of quan- 



tized high-energy vibrations in steady state fluorescence. 

Absorption and fluorescence spectra are computed in 
the absence and presence of spatial correlations and av- 
eraged over 10000 realisations of Gaussian static disor- 
der with FWHM= 500 cm~^. These results are pre- 
sented in Figure [3] together with the experimental data 
of Novoderezhkinin et a/P^. The polaron theory gives 
the correct location of the absorption peak frequency, 
though it does tend to underestimate the absorption just 
beyond this maximum (Figure Isk). The differences in 
experimental and calculated absorption spectra may in- 
dicate that full quantum treatment of selected vibrations 
may be needed. Calculated and experimental fluores- 
cence spectra are presented in Figure [SJd. We recover 
the measured experimental peak fluorescence after shift- 
ing the theoretical result by 6nm to the blue. Upon this 
correction, there is good qualitative agreement between 
theory and experiment. 

Importantly, the main features of fluorescence are only 
recovered when quantized vibrations are included. While 
the continuos component of the spectral density Jcm{^) 
determines the width of the fluorescence proflle, quan- 
tized high-energy vibrations described through Jqm{^) 
give rise to the asymmetry observed towards the blue 
(around 630 nm). The linear spectra are found to be un- 
changed as the degree of broadening, 6, of the modes is 
varied from 1 to 50 cm~^, being in good agreement with 
the results of Novoderezhkin et al?'^ . 

Finally, it is worth noting that we obtain near identi- 
cal spectra in the absence and presence of spatial corre- 
lations. The PE545 complex is a system whose excitonic 
states remain highly localized for a wide range of spatial 
correlations. Consequently, the dipole-dipole correlation 
function and the resulting spectra are insensitive to the 
degree of the spatial correlations. In contrast to the linear 
spectra, excited state dynamics is susceptible to spatially 
correlated fluctuations as we discuss in the next section. 



B. Signatures of quantized vibrations in dynamics 

In this section we highlight the effects of quantized vi- 
brations by systematically comparing dynamics with and 
without such vibrations as well as comparisons when cor- 
related fluctuations are included. We shall illustrate how 
the presence of quantized vibrations is evidenced in trans- 
fer times and pathways of energy transport in PE545 
and in the beating proflle of populations and coherences. 
Before proceeding to the dynamics, however, let us flrst 
briefly discuss two important aspects regarding the initial 
state and the basis chosen to describe excitation dynam- 
ics. 

Initial state.- We consider excitation by an ultra- 
fast laser pulse in a time scale much smaller than the 
timescale of the solvation dynamics. Under these condi- 
tions, the excitonic system will undergo a vertical Franck- 
Condon transition where the vibrational degrees of free- 
dom remains in thermal equilibrium during the fast ex- 



citation process. Therefore we may assume that the ex- 
citation pulse can excite an excitonic state |e^) of the 
bare Hamiltonian (See Eq.Q). In this work we assume 
that the laser pulse is shifted well into the blue end of the 
spectrum, such that it only excites the highest excitonic 
eigenstate of the \e^) corresponding to a highly local- 
ized excitation on PEBsqd- This initial state is chosen 
in order to investigate the transfer pathways within the 
PE54 5 com plex and to allow comparisons with previous 
worka^^I^. Furthermore this situation will let us eval- 
uate whether coherent superpositions of excitonic states 
are evolved during the dynamics even if they were not 
present in the initial state. 

Basis of Relaxation.- It is worth taking a moment to 
consider the correct basis in which to describe the re- 
laxation dynamics. Assuming the single-excitation sys- 
tem relaxes towards a quasi-equilibrium state, a natu- 
ral basis of relaxation will be that in which the den- 
sity matrix becomes diagonal as the steady state is ap- 
proached. For arbitrary system-bath coupling, such a 
basis may be constructed in terms of the eigenstates of 
an effective Hamiltonian describing equilibration^^. In 
practice, however, the evaluation of this effective Hamil- 
tonian is not a straight forward task and therefore we 
must choose an approximate basis. Our calculations in- 
dicate that the long-time density matrix in the polaron 
frame is diagonal in the basis of the renormalized Hamil- 
tonian Eq.([5|. We therefore choose such a basis to de- 
scribe relaxation but now in the untransformed frame. In 
general, these renormalized states depend on the charac- 
teristics of system-bath interaction such that it will be 
different in the absence or presence of correlations. In 
the absence of solvent-induced correlations, we flnd that 
fluctuations induced by Jcm{^) lead to all bath-induced 
renormalisations ^m^ evaluating to zero. As a result, the 
master equation EqJ7| becomes a set of non-Markovian 
rate equations describing incoherent excitation transfer 
between sites. In the presence of spatial correlations cor- 
responding to a speed of sound of 1000 ms~^, renormal- 
isations are no longer zero and quantum coherent con- 
tributions to dynamics can be observed. The basis of 
relaxation in this case is a renormalized excitonic basis. 



1. Transfer times and pathways 

Quantized vibrations drastically change the transfer 
times and pathways for excitation transfer in PE545. To 
illustrate this we begin by presenting in Figure |4] snap- 
shots of the distribution of the excitation in real space at 
different times in the presence of these vibrations (Fig- 
ure 



4a) and in their absence (Figure 4b). The main fea- 



ture that can be appreciated from these snapshots is that 
in both situations excitation energy initially localized on 
PEB50D is flnally distributed to the three lowest lying 
chromophores PEB82C, DBVa and DBVb- However, the 
timescale of this process is 5 ps when assisted by quan- 
tized vibrations and becomes 10 ps if they are neglected. 
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(a) Spatial distribution of excitation energy in presence of quantized vibrations. 
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(b) Spatial distribution of excitation in absence of quantized vibrations. 

FIG. 4. Snapshots of population dynamics. The filling of each bar is proportional to the probability of each site to be excited, 
with a maximum height denoting a population of 0.5 or above. 
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FIG. 5. Comparison of site population dynamics with and without quantized vibrations. To make clear the differences 
between the two situations only two lines have been labeled. Dashed lines represent dynamics when solvent-induced correlated 
fluctuations are considered. 



Importantly, the timescale observed in experiments is 2. Population dynamics 

about 5 ps agreeing with our predictions. 



Comparison of the snapshots at 0.1 ps indicates that 
quantized vibrations activate transfer to its nearest 
neighbour chromophore PEB50C (also the most strongly 
interacting site). This can also be appreciated in Fig- 
ure [5^ depicting site population dynamics as a function 
of time. In contrast, when no quantized vibrations are 
included, excitation is distributed to the chromophore 
nearest in energy PEB82D as can be seen in Figure [sJd) . 
One striking feature of the distribution profile at 0.7 ps in 
Figure [4a| is that excitation energy has been widely dis- 
tributed over all the complex, while in the absence of fast 
vibrations energy is concentrated mainly on protein sub- 
unit D. Again, results at these intermediate timescales 
are in good agreement with experimental observations 
irP^. The features observed in the presence of fast vi- 
brations allow PE545 to perform its biological function: 
distribution of excitation across the antenna implies that 
energy can be efficiently transferred to other nearby an- 
tennae, and fast transfer to units A and B, which are 
those probably in contact with PSI or PSlP^, will aid 
energy conversion. This is further discussed in section 
\V\ Finally, it is worth nothing that the participation 
of quantized vibrations in excitation dynamics guaran- 
tees that the general timescales of transfer and overall 
distribution of excitation energy remains practically un- 
changed as correlated fluctuations are considered (Figure 
Isk). This is no longer the case when fast modes are ex- 
cluded. In this situation, since transport is mostly driven 
by thermal fluctuations, a small degree of correlations 
slows transfer (Figure pp) . The interplay between corre- 
lated fluctuations and quantized vibrations does exhibit 
interesting features when contributions of each individ- 
ual mode to site populations are analyzed as shown in 
the next section. 



High-energy vibrations quasi-resonant with excitonic 
transitions manifest themselves not only in the activa- 
tion of transport between widely energetically separate 
states but also in the oscillatory pattern of the dynam- 
ics. This arises primarily from the non-equilibrium dy- 
namics of the vibrational components. To illustrate these 
features, we focus now on the population evolution of of 
the renormalized state 1^4), which is mostly localized on 
the central pigment PEBsoc- Figure 6jleft) presents the 
population of PEB50C when the 14 vibrational modes are 
included and also when only the modes with frequencies 
^QM = 1111 cm~^ or ujQM = 1450 cm~^ are present. 
While ujQM = 1111 cm~^ is in close resonance with the 
energy detuning between the initially excited site and the 
central site, the mode with ujqm = 1450 cm~^ has the 
largest coupling to each site. We find that both modes 
are able to activate transport in comparison to the situa- 
tion when no mode is included (curve labeled (I)), with a 
larger rate enhancement achieved with the quasi-resonant 
mode as expected. The Fourier spectrum of the popula- 
tion dynamics with the full set of 14 modes is somewhat 
complicated (not shown). Hence we investigate the os- 
cillatory behaviour when only one mode is present (see 
(I) to (HI) in Figure p[ left)). These figures show that 
the subtle oscillations of these populations are due to the 
corresponding vibrational modes. 

It is important to point out that off-resonance modes 
such as ujQM = 1450 cm~^ are able to activate transport 
to PEB50C due to the interplay with thermal-induced 
fluctuations induced by Jcm{^)' fluctuations in exci- 
tonic energy splittings cause excitonic transitions to be 
close to an otherwise off-resonance vibration. Equally, 
fluctuations bring out of resonance transitions which are 
otherwise quasi-resonant with a fast mode. The inter- 
play between thermal fluctuations and quantized vibra- 
tions bring us to the discussion of the interesting be- 
haviour observed when a small degree of correlated flue- 
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FIG. 6. Population dynamics of the excitonic state localized on PEB50C in the absence of spatial correlations (left) and in 
the presence of solvent-induced correlations (right). Lines labeled with (IV) denote population when all modes are included. 
Panels (I) to (III) depict Fourier transforms of populations when no LM is included, or when only ujqm — llllcm"^ or 
OJQM — 1450cm~^ are present. 



tuations is included (see Figure loFright)). In this case, 
excitonic transitions remain quasi-resonant with a par- 
ticular mode that fully enhances transport as shown in 
Figure |6|right) (curve labeled (II)). Furthermore, the 
population dynamics exhibits frequency components of 
both excitonic and vibrational origin. This can be seen 
in Figures loFright) (I) -(III) which show the Fourier spec- 
trum of populations when there is a small degree of 
solvent-induced correlations. When only the mode with 
^QM = 1450 cm~^ is included the spectrum exhibits 
three peaks. One peak at ~ 1450 cm~^ due to the local- 
ized vibration assisting transport and a low-amplitude 
peak at ^ 1040 cm~^ which compares very well with 



the energy difference between eigensates \as) and 1^4). 
Hence we may attribute this frequency component to a 
quantum-coherent excitonic oscilation and indicates that, 
in the basis chosen to describe the dynamics, popula- 
tions can witness coherences. The feature at 400 cm~^ in 
Figures [gJIII) (right) does not correspond exactly to any 
transition frequency, therefore we believe that this peak 
arises from a beating between a frequency of excitonic 
coherence ( uj ~ 1040 cm~^) and the vibration-induced 
oscillation of frequency uqm = 1450 cm~^. For compar- 
ison, we present the Fourier spectrum when no mode is 
present (see Figures loFright) (I)) and when only the mode 
^QM = 1111 cm~-^ is included (Figurespfll) (right)). The 
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tions are included and their frequency spectrums 



small scaling factor of the spectra in Figure JGJ^right) (I) is 
consistent with the very small amplitude of the popula- 
tion oscillations observed when vibrational modes are not 
present. 



3. Coherence dynamics 

As discussed in the previous section, quantum- coherent 
contributions to dynamics appear when solvent-induced 
correlated fluctuations are included. In this situation, 
bath-induced renormalizations do not vanish and they 
range between /328 ~ 0.11 and ^4^ ^ 10~^ such that 
quantum coherent dynamics associated to H^i is at play. 
We therefore now investigate how non-equilibrium vibra- 
tions manifest themselves in the oscillatory profile of ex- 
citonic coherences. For simplicity, we consider the situa- 
tion where only the mode ojqm = 1111 cm~^ is present 
and show in Fig. 7a the evolution of the coherences be- 
tween the two highest energy renomalised eigenstates i.e. 

Pa8,a7 = (^8 |p(^) I^^t) aS Wcll aS Pa8,a6 = (o^S |p(^) |<^6) • 

These off-diagonal matrix elements correspond to coher- 
ences investigated in room-temperature two-dimensional 
spectroscop}J^. 

Large-amplitude oscillations are observed in the dy- 
namics of both pa8,a7 and pa8,aQ but they last only for 
about 50-100 fs. On longer timescales, there is also a 
small-amplitude oscillatory component -a hundred times 
smaller than the oscillations in the short time scale. To 
investigate the origin of the two different oscillatory fea- 
tures, we once again calculate their Fourier transforms. 
The frequency spectrums over the 1 ps timescale are pre- 
sented in Figure |7b| Two features should be noticed. 
First, there is a very broad component due to the fast de- 
caying oscillations seen at short times. For the matrix ele- 
ment Pa8,a7 this broad feature occurs at 500 cm~^, which 
corresponds to the energy difference between eigenstates 
\a^) and \aj). Similarly for pa8,aQ this broad features is 
observed around uj = 600 cm~^, which corresponds to 
the energy difference between eigenstates |a8) and \a7). 
Therefore, we may attribute these frequencies to quan- 
tum coherent evolution associated to H^i. In contrast, 
the second feature within the Fourier spectra is much 
subtler and occurs in both cases at roughly 1100 cm~^. 

To investigate the fine structure of the spectrum of 
long-lived oscillations we perform the Fourier transform 
only over the time interval 250 fs to 1 ps (Figure 7c ). 
Both off-diagonal elements exhibit a peak at 1100 cm~^, 
which can be attributed to vibration-induced oscillations 
due to the localized mode with ujqm = 1111 cm~^. Inter- 
estingly, however, the matrix element pa8,aQ exhibits an 
additional peak at 1040 cm~^, an energy scale compara- 
ble to the transition frequency between eigenstates \a^) 
and 1^4). This is an indication of coherence-coherence 
coupling during the evolution and shows that excitonic 
coherences are generated during the dynamics even if 
they were not initially created. 

In summary, we have shown that in the presence of 
a small degree of correlations supporting coherent evo- 
lution of electronic excitations, the beating pattern of 
excitonic coherences in the presence of vibrations out of 
equilibrium has a dual origin: short-time oscillations last- 
ing for about 50-100 fs reflecting the time scale associ- 
ated to the electronic Hamiltonian and a superimposed 



12 



0.5 


c^-^ 


^ 




- 


^0.4 




^-- 




- 




^^~~~~~~~~-=— =r— ""^ ' " " 




CO 


/ 


B^ 


^^^ 




|0.3 




lo.2 


-""a 




w 


".-^ 


^^"^ 






0.1 


/ 


"~~^ 


■—-—____ 









12 3 4 5 

Time (ps) 

(a) Non-cascaded energy transfer between protein subunits when 
quantized vibrations are included 
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(b) Cascaded energy transfer from protein subunit D to A when 
quantized vibrations are neglected 

FIG. 8. Total population within each protein subunit of 
PE545. For subunit D we have neglected the contribution 
of the initially excited chromophore PEBsod- 
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FIG. 9. Performance measures in the presence (I) and absence 
(II) of quantized vibrations. Insets illustrate the behaviour 
when solvent correlations are considered. The dotted line at 
1 ps indicates the estimated transfer time between antennae 




vibration-induced modulation that survives close to the 
picosecond timescale. Importantly, the two timescales 
identified here compare very well with recent experimen- 
tal results on a similar complex^. 



V. ROLE OF QUANTIZED VIBRATIONS 

The results of the preceding section allow us to clearly 
point out three key roles of the strong-coupling between 
highly localized electronic excitations and quantized vi- 
brations in the biological function of PE545. 



A. Promotion of non-cascaded energy transfer 

In the absence of quantized vibrations exciton-phonon 
interactions are dominated by the continuous distribu- 
tion of modes that induce transitions between states near 
in energy. Given the local nature of excitonic states, 
transfer dynamics follows the hierarchy of energies of 
the chromophores. In light of the intimate connection 



between the hierarchy of energies and the physical ar- 
rangement of the pigments in the protein subunits (Fig- 
ure Ic), we conclude that when quantized vibrations are 



neglected, electronic excitation moves through PE545 in 
a sequential manner, from D onto subunit C and finally 
to subunits B and A as shown in Figure [Sb] Clearly this 
would be a very unfavourable dynamics for the biological 
function of transferring energy to other antennae. 



In contrast, quantized vibrations take advantage of 
electronic interactions to promote non-sequential trans- 
fer by bridging the energy gap between subunits D and 
C and spreading out energy from each of them to A and 
B as shown in Figure |8a| In short, quantized vibrations 
in close resonance with excitonic transitions provide a 
built-in mechanism to control energy transfer pathways. 
In general, explotation of environment and system reso- 
nances is an important design principle for directing en- 
ergy transfer^^J 
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B. Optimal spatial and temporal spread of excitation 

High-energy vibrations in PE545 guarantee an effective 
distribution of energy throughout the complex towards 
the acceptor states. The activation of energy transport 
between PEB50D and PEB50C generates both wider spa- 
tial spread of excitation to allow transfer between anten- 
nae and also faster transfer to the final acceptor states 
(localized on DVB a, DVBb and PEB82c) to guarantee 
energy conversion if an PE545 antennae is in close prox- 
imity to complexes like the PSII or PSI. 

Quantitative support of the above statement is shown 
in Figure [9] which summarizes the dual role of the quan- 
tized vibrations in promoting faster transfer and wider 
spatial distribution of excitation energy. The total popu- 
lation of all the acceptor states is depicted in Figure [9a 
and clearly shows a significant enhancement of the trans- 
fer rate in the presence of promoting vibrations. To quan- 
tify the spatial spreading of excitation energy in PE545 
complex we compute the variance cr^(t) in the spatial 
position of the excitation defined as follows: 



(j\t) = ^Pm{t){\r, 



\t)). 



(17) 



where Pm{t) denotes the population on chromophore m at 
a given time t, r^ is the position vector of chromophore 
m and r{t) is the mean position of the excitation. The 
dynamics of this measure of spreading is shown Figure 
9b Using generalized Forster theory (details of which are 
not presented here) we have estimated that transfer times 
between adjacent PE545 complexes is of the order of 0.7 
to 1 ps. In this biological timescale, quantized vibrations 
guarantee efficient inter-complex transfer by promoting 
a wider distribution of excitation energy (dotted vertical 
line in Figure 9b) and faster transfer to acceptor states 
(see dotted vertical line Figure 9a). These conclusions 
also hold when a small degree of correlated ffuctuations 



is included as shown in insets of Figures 9a and 9b 



C. Exploitation of quantum coherent dynamics to tune 
useful resonances 



We have discussed how a small degree of correlated 
ffuctuations supports quantum-coherent dynamics of ex- 
citonic states and ensures tuning of a resonance with a 
preferred vibrational mode that enhances energy trans- 
port to a targeted site in the complex. For clarity this 



is summarized in Figure IQa showing the enhancement 
of the population of the central pigment PEB50C when 
coherent dynamics (correlations) is turned on and only 
the mode ojqm = llllcm"^ is considered. These results 
lead us to argue, that in the presence of quantized vibra- 
tions, any mechanism supporting coherent evolution of 
excitonic states, will also ensure that resonances enhanc- 
ing transport are activated. The rationale behind this 
idea is simple: coherent evolution effectively implies that 



excitonic gaps during dynamics do not ffuctuate much 
and hence frequencies associated to coherent evolution 
can be synchronized with quasi-resonant vibrations for 
longer times, which in turn enhances population trans- 
port. 

To illustrate this idea further, we consider now the sit- 
uation where quantum-coherent mixing of electronic and 
vibrational degrees of freedom is allowed, which generates 
room-temperature coherent evolution of excitonic states. 
To show this we consider the full Hamiltonian evolution 
of an enlarged system that inlcudes both electronic states 
and quantized vibrations. For simplicity we consider a 
simpliffed model consisting of the central dimer formed 
by PEB50C and PEB50D with each pigment coupled to 
only one mode of frequency ujqm = llllcm"^. As men- 
tioned before, this mode is in close resonance with the 
energy difference between the two excitonic states of the 
dimer which we denote X and Y. To capture the effects 
of the external thermal background we consider each site 
to be subject to a small amount of pure dephasing. De- 
tails of this simpliffed model are presented in Appendix 
[B] The initial state of the dimer-plus-modes system is 
separable, with the electronic part in the highest energy 
exciton state X (and hence no initial exciton coherence) 
and the high-frequency modes in thermal equilibrium. In 
this way, if no modes are present, there is no transport. 
In this situation we see that coherent exciton-mode inter- 
actions: (i) activate population transport to the low-lying 
excitonic state Y and (ii) create and modulate long-lived 
(picosecond) excitonic coherences (Figure 10b). Creation 
of excitonic coherences is, of course, expected as the 
initial state is not an eigenstate of the full dimer-plus- 
modes system. Notice that our calculations are at room- 
temperature where the thermal population of the mode is 
mostly limited to the lowest number states, hence the im- 
portance of its quantum behavior. The full advantages of 
quantum features of high-energy modes in transport pro- 
cesses will be soon presented in a separate publicatiorP^. 
A recent worlP^ has also considered modulation of exci- 
tonic coherences in the Fenna-Matthews- Olson complex 
due a vibration whose energy is comparable to the ther- 
mal energy scale (180 cm~^). It would be interesting to 
compare the differences that can arise in exciton coher- 
ence dynamics and non-equilibrium transport properties 
depending of the energy scale of the vibration involved. 



VI. CONCLUDING REMARKS 

The role of vibrations for adjusting and enhancing 
transport in wide variety of biologic al and chem ical sys- 
tems is currently of central intereslP'^^^^^^'^^^. In this 
article we have illustrated the importance of quantized 
vibrations for efficient energy transport in photosynthetic 
systems where excitonic states are highly localized such 
as the PE545 antennae protein present in cryptophyte 
algae. The main feature allowing these systems to take 
full advantage of fast vibrations is the near-resonance be- 
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(a) Enhanced population of PEB50C when coherent dynamics 

due to correlations is present (dashed line). Solid line indicates 

non-Markovian yet incoherent dynamics. In both cases only 

mode UJQM = llllcm"-'^ is considered 



antennae. We are currently pursuing such study and 
our preliminary results confirm the features presented in 
this article and indicate that rich nonclassical and non- 
equilibrium behavior arise in these systems^^. Atomistic 
calculations identifying the structural origin of such vi- 
brations are therefore crucial. From the experimental 
view point it is important to develop more direct mea- 
surements of the influence of such vibrations in chan- 
nelling and/or enhancing excitation energy transfei^ as 
well as conclusive signatures of their quantum properties 
in conditions of biological relevance. The integration of 
physical descriptions as the ones presented in this arti- 
cle, together with atomistic calculations and experimen- 
tal characterization will provide useful insights for the 
design of antennae prototypes that exploit electronic in- 
teractions and promoting vibrations for controlled and 
enhanced performance. 
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FIG. 10. Fast vibrations allow exploitation of quantum co- 
herent evolution of excitonic coherences to enhance transport 
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Appendix A: Linear Spectra Within the Polaron Formalism 

In this Appendix, we derive an expression for the 
dipole correlation function within the polaron formalism 



used in section [III| The dipole-dipole correlation function 
is defined as 



tween vibrational frequencies and the energy gaps be- 
tween excitonic states. An important implication of this 
close energy match is that quantum-coherent evolutions 
can then tune the resonances and direct excitation en- 
ergy in a preferential manner. Therefore, we argue that 
in photo synthetic complexes where transport is vibration- 
activated, a fundamental biological function for quantum 
coherent contributions to dynamics is to support reso- 
nances that promote fast and effective energy distribu- 
tion. 

The non-Markovian master equation in the polaron 
framework captures some important effects of the non- 
equilibrium dynamics associated to the vibrations of in- 
terest such as long-time modulation of the excitonic co- 
herences at room temperature, predicting time scales 
in good agreement with those reported in experiments. 
However, a full quantum treatment of the interactions 
between excitons and these quasi-resonant vibrations is 
needed both to assert the accuracy of the results here 
presented and to have a more complete picture of the 
rage of non-equilibrium effects in different chryptophyte 



Cd-d{t) = Tvs+BMt)KO)pB\0){0\} (Al) 

where /i(t) = e~*^^ /i e*^^ and the system dipole operator 
/i is defined as the sum of the dipole operators for the 
individual chromophores as follows: 



/i = };^/i^|0)(m|+/i;,|m)(0|. 



(A2) 



Let us now transform into the polaron frame de- 
fined by the transformation H = e^He~^ ^ where S = 

Em ^m^^ Ek(^k,m^k " ^k,m^k), with /lk,m = ^k,mM. 

Within the polaron frame, the dipole-dipole correlation 
function becomes: 



Cd-d{t) = Trs+Bmt)mpB\0){0\} (A3) 

while the dipole operator in the polaron frame is given 
by: 
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fi = Y, (n^("V™)Mm|0)(m| + /i.e.) 

m k 

= XI ( n Di-Km)l^mUma\0) {a\ + h.c}j , (A4) 



cx^m k 



where D{h\^j) = e^^'^^k-^k,i^k ^g ^-j^^ ]^^-^i^ displacement 
operator of mode k due to interaction with site j. 

Let us now proceed to calculate the dipole-dipole cor- 



relation function in steps. Firstly consider: 

fl{0)pB\0){0\ = ^ (n^(-^k,m)Mm^ma|0)(a| + h.c}j\0) {0\pB 
a,m k 

= ^n^(^k,m)Mm^ma|<^)(0|pB. (A5) 



a,m k 



Then the dipole-dipole correlation function is: 



Cd-d{t) =^^]JTrs^B{e-''''[D{-h^^r^)PnUn^\0){P\^D{h^^r^)^>^^ 
cx.m /3,n k 

o^.m /3,n k 

X e'**Z)(/ik,m)Ai>ma|a)(0|pBe-^**} 
ZZ n Trs+B{(D{-h^,,^)^i^Unp\Q) (/?| + Z)(/ik,„XM„/3 1/?) (0|) 



oc,m (3,n k 



xe^^^I^(V-)e"^^'e^^V;.^m«|c.)(0|pBe-^^^}. 



(A6) 



I 

At this point, we make a number of approximations where we have defined the dressed bath correlation func- 
to simplify the expression for the dipole-dipole correla- tion Ca{t) as 
tion function. Firstly, we shall assume that the interac- 
tion Hamiltonian Hj does not perturb the bath. This 

amounts to assuming that the bath is in the stationary, Ca{t) = y^Pml-^nUmaUnaTrB{D{—h\^^ri)D{h\^^rn{t))PB}- 
thermal state within the polaron frame {pb) at all times. mn 

Under this assumption we may write 



(AlO) 

To complete the calculation of the linear spectra we 

need to know how each of the off-diagonal element p^o 

e'^^D{h\^^rn)e~'^^ = e*^o^i:>(/ik,m)e"*^°^ = D{h\^^rn{t))- evolves with time. To do so, let us consider the Marko- 

(A7) ^^^^ secular polaron master equation in the interaction 
picture: 



The second assumption we shall make is that all system 
coherence terms decouple. This amounts to making the 
secular approximation for the system dynamics. There- 
fore, we may write: 



e''''^'^'pl,UrnaW){{)\pBe-'^'^'^' = pl,UrnaPa0(t). (A8) 

Under these assumptions, the dipole-dipole correlation 
function has the following very simple form: 



Cd-d(t) = ^ Trs{p*^PnUmaUn^\^){P\pao{t)} 

a(3,mn 



dp{t) 
dt 



XC«(t)Pao(t). 



(A9) 



+r(s,^, [5^., Sim] -^Ti%^ [si, Sim] 

+ h.c.). (All) 



Details of this master equation can be found in Kolli et 
alP^. The equation of motion for the off-diagonal element 
Pao is then given by: 



PaOyt) — 7 ^ ( -L afi,fia~^^ /ia/ia^-'- afi,aiJ.~^^ iJ.a,afi j PaO[t) • 

(A12) 

The general solution for the evolution of the coherence 

terms, in the Schrodinger picture, is then simply given by 
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Pao{t) = exp(icc;c^t)exp(— rQ,t)pc^o(0)7 where have defined 



the total dephasing rate as Fq, = ^ ( ^a/L,iJ,a + T}^ 



r(3) 



^(4) 



. Therefore, the final expression for the 



dipole-dipole correlation function becomes: 



Cd-d{t) = Y,Ca{t)exp{-t^t)exp{iuJat). (A13) 



Appendix B: Dimer model 

This model considers the central dimer in PE545 
formed by pigments PEB50D and PEBsoc- The exci- 
tonic Hamiltonian is given by H^ = eia^a^ + 620-2 0-2 + 
V{aia2 + cr2"crf ) where the strength of the electronic 
interaction is F = 92 cm~^ and the exciton eigenstates 
are labeled X and Y. Each pigment is coupled to a lo- 
calized mode with frequency ujqm = HH cm~^. The 
free Hamiltonian for the modes is i^modes = ^qm (^1^1 + 
62^2) while their linear interaction with the dipole- 
coupled dimer is given by J^g-modes = 90'iO'i{bl -\-bi) -\- 
QG 2 (J 2(^2 + ^2)- Here the strength of the site- mode cou- 
pling is ^ = 267 cm~^. The initial state of the full dimer 
plus modes system ^(0) = /)e(0) (8) Pmodes(O) is propa- 
gated under the total Hamiltonian H = H^ -\- i^modes + 

-^e— modes- 

The initial state of the dimer is chosen to be the 
highest energy eigenstate X such that in the absence of 
the harmonic modes no transport occurs. The modes 
are assumed to be initially in thermal equilibrium at 
room temperature (300K). For the purposes of calcu- 
lations, the Fock space of the modes is truncated to 
a maximum occupation M = 8 which fully captures 
the dynamics of the high-energy modes for the param- 
eters considered. We also introduce pure dephasing 
of onsite energies which results in a dissipator of the 
form D{g) = ^Yjii^t^iQ^t^i - \Wt^i^Q}) with 
7 = 50 cm~^. 

See supplementary material for details on the 
polaron formalism and a table of renormaliza- 
tion factors of electronic couplings when solvent- 
induced correlations have been included 
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